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1. Introduction 



Abstract. We establish the Lagrangian nature of the discrete isospectral and isomonodromic dynamical 
systems corresponding to the re-factorization transformations of the rational matrix functions on the Rie- 
mann sphere. Specifically, in the isospectral case we generalize the Moser-Veselov approach to integrability 
of discrete systems via the re-factorization of matrix polynomials to a more general class of matrix rational 

functions that have a simple divisor and, in the quadratic case, explicitly write the Lagrangian function for 
such systems. Next wc show that if we let certain parameters in this Lagrangian to be time-dependent, 
the resulting Euler-Lagrangc equations describe the isomonodromic transformations of systems of linear 
difference equations. It is known that in some special cases such equations reduce to the difference Painleve 

^ equation. As an example, we show how to obtain the difference Painleve V equation in this way, and hence 

O we establish that this equation can be written in the Lagrangian form. 

in 

The theory of completely integrable systems and soliton equations is justly known for its rich and often 
^ r*^ unexpected connections with a wide range of other branches of mathematics and mathematical physics. 

In recent years its discrete variant, the theory of discrete completely int{\Q;rabl(; systems started to attract 
a considerable amount of attention. This subject is a part of a more general field of discrete Lagrangian 
I I mechanics, which itself is getting more and more important partly due to the development of new numerical 

algorithms based on discrete variational integrals, see, for example, the recent survey by J. Marsden and 
M. West, [MWOl]. In a series of papers [Ves88, MV91, Ves91] A. Veselov and J. Moser showed that the 
discrete analogues of many classical integrable systems, e.g., the Neumann system and the spinning top, are 
related to the re- factorization transformation of certain matrix polynomials. Such a representation explains 
V/~j the integrability mechanism for these systems, since it is a discrete version of the Lax-pair representation 

and, similarly to the continuous case, it can be used to integrate the system using theta functions. A large 
number of such examples can also be found in a recent encyclopc^dic book by Suris, [Sur03]. The relationship 
T-H between matrix factorizations and integrable systems was observed earlier by Symes [Sym82], see also a 

1^ related work by Deift et al, [DLT89]. 

In this paper we construct a generalization of the Moser-Veselov approach from the matrix polynomials 
^ to a general class of rational matrix functions whose determinant divisor is simple. On one hand, such a gen- 

eralization is important if one wants to consider discrete integrable systems that have higher-genus spectral 
curves. On the other hand, these rational matrices play an important role in the theory of isomonodromic 
transformations of systems of linear difference equations recently developed by A. Borodin [Bor04] , see also 
[Kri04]. Similarly to the continuous case, under certain conditions such transformations give rise to the 
discrete Painleve equations from Sakai's hierarchy, [SakOl], thus clarifying the geometry of these equations. 

Our main result is that both the isospectral and the isomonodromic dynamics generated by the re- 
factorization transformations on the space of such rational matrices are Lagrangian. We focus our attention 
on the two-pole case and, by introducing a special coordinate system, explicitly write down the expression 
for the Lagrangian functions. We also show that in the isomonodromic case the resulting dynamics for the 
rank two rational matrix functions with two simple poles gives rise to the difference Painleve equation dPV 
of the Sakai's hierarchy, thus establishing that this equation can be written in the Lagrangian form. 
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In the remainder of the introduction we give a detailed description of the setup of the problem and of our 

results. 



1.1. Discrete Lagrangian Systems. Continuous dynamical systems can be considered in the Lagrangian 
or in the Hamiltonian framework. Of those two, the Lagrangian approach is the one that naturally generalizes 
to the discrete case, see [Ves88, MV91, MWOl]. Let Q be the configuration space of our system and let 
n e Z be the discrete time parameter. In the continuous case the Lagrangian £ € T{TQ) is a function on 
the tangent bundle of Q. For the discrete case we need to change the point (Q, Q) in the tangent space to 
the pair of points (Q, Q) in the configuration space itself . The Lagrangian then becomes a function on the 
square of the configuration space, C €: T{Q x Q). The action functional <S is then defined on the space of 
sequences {Qfe}, fc € Z by the formal sum 

'5({Qfe}) = I]'C(Qfe,Qfe+i), (1.1) 

k 

and the variational principle 5S = Q that selects the trajectories of the system, when written in a coordinate 
chart, takes the form of the discrete Euler- Lagrange equations 

|^(Q,Q) + ||(Q,Q) = o, (1.2) 

where we use the notation Q = Qfe_i, Q = Qfc, and Q = Qfc+i. These equations then implicitly define the 
map (or, more precisely, a correspondence) Q = ^(Q,Q), which, in turn, defines the shift (or step) map 

^:QxQ^QxSby <?(Q, Q) = (Q, Q). The map is symplectic w.r.t. the 2-form a = — — dX A dY 

on Q X Q. Alternatively, defining the discrete version of the Lcgcndrc transform by defining the conjugated 
dC 

momentum P = ^^(Q) Q) ^ IqQ, we see that the discrete Euler-Lagrange equations are equivalent to the 

system 

P = -^(Q,Q) 
P = ^(Q,Q), 

where the first equation follows from (1.2) and is an implicit equation for Q = Q(Q,P). Hence we get a 
map ip' : Tq T~, !f'(Q,P) = (Q,P) which is symplectic w.r.t. the standard symplectic structure. In what 
follows by the equations of motion of a discrete Lagrangian system we mean either (1.2) or (1.3), with the 
corresponding discrete dynamics given by the maps or tf' respectively. 

1.2. Discrete Integrable Systems. The Moser-Veselov approach to the discrete integrable systems is 

based on the discrete version of the Lax pair representation and can be briefly described as follows. Given 
a discrete dynamical system, we look for a class P of matrix polynomials L(z) in a spectral variable z and 
a parameterization map r] : Qx Q ^V, defined on some dense open set, such that: 

(i) there is a well-defined factorization rule ^{z) = L,i{z)L,2{z) , where the ordering is important, such that 
the re-factorized matrix L obtained by the interchanging the order of the factors, L(z) = L2(2;)Li(z), 
is again in V and hence can be written as L(2:) = Li (2)1:2(2:), this rule defines the re-factorization map 
R:V-*V] 

(ii) under the parameterization r] the re-factorization map R corresponds to the shift map $ of our discrete 
dynamical system. 

Note that the re-factorization map can also be written in the form 

L = M-^LM, (1.4) 

where M = lii{z). Equation (1.4) is known as a discrete Lax pair representation of the system. Similarly 
to the continuous case, finding such a representation shows that the dynamic of the system is isospectral. 
Therefore it preserves the spectral curve F of the operator L(z), which implies the integrability of the system 
and also makes it possible to obtain the ^-function formulas for solutions of the system in the usual way, see 
[Ves91]. 



1.3. Rational Matrices Anzats. Let us now allow L(z) to bo a nicroniorpliic; r x r-matrix function on 
the Riemann sphere with the poles zi, . . . ,Zn- We restrict our attention to the matrices that are generic in 
the following sense: 

(i) all poles Zf of L(2) are simple; 

(ii) the divisor of L(z) is simple as well, where by the divisor o/L(z) wc mean the divisor of its determinant 
function, T> = (L{z)) = (det L(z)) = J2i ~ J2j Q- This is equivalent to the condition that the residue 
matrices Lj, := res^j^ L(2;) are of rank one. 

For a fixed divisor 2? wc denote the space of all such matrices by A4^ . 

Without any loss of generality we can further restrict out attention to the case zq = oo ^ V and Lq = 
lim^^oo L(2) is invertible and diagonalizable. Any such matrix L(2;) has two different representations, 
additive: 

L(^) = Lo + ^^^, (1.5) 

^ Z-Zk 



and multiplicative: 



k 



Zk 



where Gfe is a matrix of rank one and A* = Lq; we are mainly interested in the multiplicative representation. 

We call the factors ^f'{z) = ^^"1 — ^^^^ multiplicative representation the elementary divisors of 

L(^). Note that the ordering of the poles determines the ordering of the factors in the multiplicative 
representation, which is what we need to define the re-factorization map. 

In this paper we restrict our attention to the two-pole case. However, since any permutation is a compo- 
sition of elementary transposition, the general case is generated by the transformations described below. 

1.4. The Re-factorization Transformation and the Isospectral Discrete Dynamical System. Let 

us now fix the divisor V = zi + Z2 — d — C2, where all four points are finite and distinct, and consider the 
re- factorization map R : A4^ —^M.^ given by 

L(^) = Bf (^)B^(^) ^ t{z) = B^(0)Bf (^) = Bf (^)B^(0). (1.7) 

Wc want to determine whether there is a natural discrete dynamical system for which this map is a discrete 
Lax pair representation, and if so, what is the Lagrangian of this system. Note that this setting is rather 
general, since no restrictions on the rank r are imposed. First it is necessary to identify a configuration space 
Q such that there is a parameterization map r/ : Q x Q — > M.^ satisfying the following diagram: 

(Q,Q) ^(Q,Q) (1-8) 



Y 

L(^) = B^(z)Bf (z) = Bf (z)B^(^) ^ L{z) = B^(^)Bf (z) = Bf (z)B^(z) . 

From this diagram it is clear that half the data in L(z) should come from Q and half should come from Q 
and moreover, this data should be of the same type to be compatible with the shift map ^. We know that 
the elementary divisors 'Rf^{z) completely determine 'L{z), and each elementary divisor is in turn determined 
by the rank-one matrix Gj = PiqJ, where p, and qj are defined up to a common scaling constant, and this 
constant can be recovered from the divisor T>. Thus, L(z) is completely determined by either pi, p2, q|, q2 

or Pi, P2, qi, q2- In view of that we take Q — (pi, q2), Q = (Pi^ql)- A priori Q € C" x (C)^, but in fact the 
resulting expressions are homogeneous in Pi, qj, and so the correct configuration space is Q = P''-! x (P''^^)'''. 
In Theorem 3.1 we give an explicit description of the parameterization map ry : Q x Q ^ ■ compute the 
corresponding equations of motion, and show that these equations of motion are the discrete Euler-Lagrange 
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equations with the Lagrangian function C, given by 



£(X, Y) = (^2 - zx) log(xtxi) + (01 - C2) log(xtA-Vi)+ 

(C2 - Ci) log(yjA-Vi) + (Ci - ^2) log(ytA-ixi). 



(1.9) 
(1.10) 



1.5. Discrete Painleve Equations and the Isomonodromic Transformations of the Systems of 
Linectr Difference Equations. Another natural discrete dynamics that can be considered on our space of 
matrices is the isomonodromic discrete dynamical system. One of the reasons this system is interesting is 
its relationship to the theory of the discrete Painleve equations. 

The recent surge of interest in the discrete version of the famous Painleve equations is in part due to the 
fact that these equations appear in the calculation of discrete gap probabilities in the theory of (determinantal) 
Random Point Processes, [Bor03, BD02]. In addition, H. Sakai in [SakOl] described a very elegant and purely 
geometric approach to the discrete Painleve equations using the Cremona action on algebraic surfaces. More 
information about the current progress in the theory of discrete Painleve equations can be found in [GR04] . 

In the continuous case there is a well-known relationship between the isomonodromic transformations of 
the flat meromorphic connections on the Riemann sphere and the Painleve equations. Thus, it is natural to 
expect that discrete Painleve equations should be related to isomonodromic deformations of matrix linear 
difference equations. However, there is a serious obstacle — the notion of monodromy for a differential 
equation has no obvious generalization to a difference equation, and only recently some significant progress 
was made in this direction. Recall that the general theory of matrix linear differential equations 



goes back the works of George Birkhofi^, [Birll]. First step in Birkhoff's approach was to use a special gauge 
transformation to clear all poles of Li{z) and make it a polynomial in z. Note that as a result we get a pole 
of higher order at infinity. Next, Birkhoff showed that there are two canonical meromorphic solutions ^i{z) 
and ^r{z) that have the prescribed asymptotic behavior for 3?(z) <C and di{z) 3> respectively. Then the 
analogue of the monodromy map is just the connection matrix C of these solutions, C{z) = '4'~^(z)\I';(z). 
Birkhoff also showed that in this situation the isomonodromic transformations L(2;) 1-^ L(2;) that preserve 
C{z) are given by 



where R{z) is a rational matrix. A. Borodin, in [Bor04], constructed a general theory of such transformations 
for polynomial L(z) and showed that it give rise to difference Schlesinger equations. These equations, when 
the dimension of the space of parameters is equal to two, can in turn be reduced to difference Painleve 
equations. In a follow-up paper [AB06], D. Arinkin and A. Borodin showed, using a more geometric language 
of rf-connections, that for some special cases Sakai's surfaces can be identified with the moduli space of 
such d-connections and the isomonodromic transformations can be though of as elementary modifications 
of d-connections, which are in turn given by difference Painleve equations (examples considered in this 
paper are dPV and dPVI). This result helps to explain the geometry behind the difference isomonodromy- 
Painleve correspondence. Recently Arinkin and Borodin found the description for the tau- function of discrete 
isomonodromy transformations for both polynomial and rational cases, see [AB07]. 

A different approach to the notion of the monodromy of a linear difference equation was suggested by 
I. Krichever in [Kri04]. In his approach the matrix L(z) belongs to the same anzats as we consider in this 
paper — it is regular at infinity, all of its poles Zi are finite and simple, and res^^ L(z) must be of rank one. For 
such matrices Krichever introduced the notion of local monodromy that can be thought of as a monodromy 
corresponding to the path around a pole, and also constructed the isomonodromy transformations, that also 
have the form (1.12). He also showed how to generalize the theory from rational to elliptic functions. 

In the present paper we consider a special case of the transformation (1.12) that has the form 



L(z) = Bf (z)B^(z) ^ Uz) = Btiz)B^iz) = B^(z + l)Bf (z) = B^(z + l)L(z)(B^(z))-\ (1.13) 



where 'Bf'{z) are the elementary divisors defined earlier. We show that, similarly to the isospectral case, 
these transformations can be written in the Lagrangian form. The main new feature of the isomonodromic 
approach is the fact that such transformation changes the divisor D to the divisor D, where Zi = Zi, (i = (i. 



*(z+l) =L(z)*(0) 



(1.11) 



L(z) = R(z+ l)L(z)R-i(z), 



(1.12) 
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Z2 = 22 — 1 , and C2 = C2 — 1 • Thus we need to make the Lagrangian C time-dependent by putting Z2 {t) = Z2 — t 
aiidC2W = C2-i: 

£{X,Y,t) = {Z2{t) - zi)log(x|xi) + (zi - C2(i))log(x^A-Vi)+ (1-14) 
(C2(t) - Ci) log(yjA-Vi) + (Ci - Z2{t)) log(y|A-ixi). (1.15) 
The time-dependent discrete Euler-Lagrange equations 

^(Qfc-i, Qk, k-l) + g^iQk, Qfe+i, k) = (1.16) 

then describe the isomonodromic dynamics (1.13). Finally, we verify, essentially following [AB06], that in the 
rank-two case equation (1.13), when written in the so-called spectral coordinates, reduces to the difference 
Painleve equation dPV of the Sakai's hierarchy, thus establishing that this equation can be written in the 
Lagrangian form. 

1.6. Organization of the Paper. In Section 2 we study properties the elementary divisors, and obtain the 

description of the re-factorization map. In Section 3 we establish the Lagrangian structure of the isospectral 
dynamics, and in Section 4 we extend this to the isomonodromic case. 

1.7. Acknowledgments. I am very grateful to L Krichever for many helpful and interesting discussions. 

2. Elementary Divisors 

In representing rational matrix functions in the multiplicative form we take each factor to be a matrix of 
the following simple type. 

Q 

Definition 2.1. An elementary divisor with the simple pole at Zi is a matrix of the form Sf-(z) = A-\ 

Z- Zi 

where Gj is a matrix of rank one and A is some fixed constant non-degenerate matrix (which is usually taken 
to be diagonal). When A = 1, we denote 'St\{z) by simply 

In this section we describe certain useful properties of elementary divisors, and also explain our normal- 
ization conventions. 

2.1. Rank-One Matrices cind Normalization. Let us first make some remarks about matrices of rank 
one. Any such matrix has the form G = pq^ for some column vector p and some row vector q^, where 
the vectors p and q^ are defined up to a common scahng constant. To explicitly keep track of such scaling 
constants during computations we need to normalize these vectors in some way. 

< Notation: Given the choice of a normalization, we denote by [v] the normalization of a vector v and by 
v{y) its normalization constant w.r.t. this normalization. Thus, v = :^(v)[v]. We also use the notation [v] 
for normalized vectors. Then any matrix of rank one can be written as G = pqt = !^(G)[p][qt] = A[p][qt], 
where A = i^{G) = i'(p)i'(q^) is the normalization constant for G. > 

For our purposes it is most convenient to work with some linear normalization. Such normalization has the 

property that any linear relation for normalized vec;tors implies the same linear relation for the coefiicients; 
if a[u] — ^k[vk], then a = J2k ^k- For example, the normalizations 

(a) E,(pr = E,(q^), = i, 

(b) (p)' = (q^)j = 1 for some choice of indexes i and j 

satisfy this requirement. 

For an elementary divisor we often use one of the following three natural normalizations: 

z — Zi y z — Zi J y z — Zi J 

where the superscript p in Af indicates that instead of normalizing p we normalize A~^p; similarly, A^^ 
corresponds to the normalization of q| A~^. 



2.2. Properties of Elementary Divisors. The following Lemma is a key technical tool for working with 
elementary divisors. 

G • 

Lemma 2.1. Let'Bf-{z) = AH — and define Q by the equation tv{GiA~^) = Zi — Q. Then the following 

Z - Zi 

holds. 

(i) dctBf (z) = ^^det A and (Bf (z))-i = A-^ (a - 

(a) Knowing how B^(z) operates on row (resp. column) vectors and also knowing the column (reap, row) 

vector of the rank one part, allows us to determine B^(z); 

• ifv = Bf (z)w, then Gi = A ( {zi - z)^ + {z - d) f^^^ ] qj; 



• «/vt = wtBf (z), then G, = p, - z)^^ + {z - Q) ^tA-ip J ^' 

Proof. To prove part (i), note that 

detBf(z)=det(l + ^^^)detA= (l + ) det A = det A, (2.2) 

since GjA"^ is a matrix of rank one. The formula for the inverse matrix can be checked by the direct 
calculation. 

To establish part (ii) we normalize the elementary divisor. Then, using the linearity property of the 
normalization, the equation v = Bf^(z)w can be written as 

[w](z-z,) + Ar[A-ip,][qt][w] 



fA-^vl = 



. .^ Af[A-V.][qt][w] 



(z-z,) + Af[qt][w: 



(2.3) 



vpr„tir,„i _ [ql]M(2-Ci) 



Multiplying both sides by [qj] gives [z — Zi) + A?[qJ][w] = — , and so 

[qT][A-iv] 

^ ' V M[qi] [qi][A-iv] ) 

Since the expression in the parentheses is homogeneous, we can remove the normalization brackets to get 
the desired result. Second formula is obtained in a similar way. □ 

2.3. The Re-Factorization TrEinsformation. Consider now the following question. Let 

B^(z)Bf(z)=Bf(z)B^(z). (2.5) 

where Bf-(z) and Bf-(z) are elementary divisors with the simple poles at Zj. What are the relationships 
between the vectors that form their rank-one paxts? To begin with, note that taking the determinant results 

in the equation \ ^^^^ = \ ^^^7 where zi = Zi by definition, and so wc must have cither 

(Z-Zj)(z-Z2) (Z-Zl)(z-Z2) 

Ci = Ci (the general case) or = ^2 and ^2 = Ci (which is a special case, since it requires a non-trivial 
relationship between the poles and the rank-one parts of the elementary divisors, 

Z2-zi= tr(G2A-^) - tr(GiA-i) = tr(G2A-^) - tr(GiA-i). (2.6) 

The theorem below explains the general case, and the special case can be considered in the exactly same 
way. 

Theorem 2.2. Let B^(z)Bf (z) = Bf(z)B^(z) and Q = Q- Then the following holds. 

(i) The vectors Pi, Pi, qj, qj are related by 

[pi] = [B^(zi)Pi] = [AB^(Ci)A-ipi] [p2] = [Bf (Z2)P2] = [ABf (C2)A-ip2] (2.7) 

[ql] = [qlB^(zi)] = [ql A-iB^(Ci)A] [4] = [q^Bf (Z2)] = [q^-^Bf (C2)A] . (2.8) 
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(ii) The vectors pi, pi, q^, q2 completely determine the elementary divisors Sf-{z), B^(^) via 

Gi = f (.1 -z,)^ + {z, - CO^I^) A, (2.9) 

V q2Pi q2A-ipi / 

= C.) + (C. - CO ^f^) A. (2.10) 

V q2A-ipi q2A-2pi J 

G2 = A f (.2 - .i)P^^ + (.1 - G)^^) , (2.11) 

V q2Pi q^A^ipi / 

Proof. The equations in part (i) are obtained by taking the residues of the equation B^{z)'Bf-{z) = 
B^(z)B^(z) at the points Zi and the residues of the inverse equation at the points Q and then equat- 
ing the normalized row and column vectors of the resulting rank-one matrices. Using Lemma 2.1 we then 

obtain the equations in part (ii). □ 

3. The Isospectral Case 

We are now in the position to describe the equations of motion for the isospectral dynamics. 

Theorem 3.1. Let Q = (pi,q2), Q = (pi,q2)> Q = (Pi)q2)> where each vector is considered modulo 
re-scaling, 

r?(Q, Q) = L(z) = B^(z)Bf (z) = Bf (z)B^(z), (3.1) 

and 

7j{Q, Q) = L(^) = B^(^)Bf (z) = Bt{z)B^{z). (3.2) 

Then 

(i) The parameterization r]{Q,Q) is given by 

,(Q,Q)=LW= (a+-2_ J,,_C,)-g4-^ + K,-C,)|^]] X (3.3) 

(A+^f(.-CO^ + (C.-<.)^)). (3.4, 

V z^Z2 \ q3>A-ipi q^A-2pi J J 

(ii) The equations of motion (Q,Q) = ^(Q,Q) have the implicit form 

( - CO^^ + (C. - C.)^^) -L- + (.. - C.)^^) (3.5) 

V q2A-ipi q;^A-2piy \^ qipi q^A-^pi^ 

f (^1 - C2)^^ + (C2 - Cl)^^) = (izi - Z.)^ + (Z2 - Cl)4rT-] ' (3.6) 

V q2A-ipi q^A-2piy y q^pi q^A-ipi^ 

and since we are only interested in the spaces spanned by pi and qj, we can take pi and q2 to be given 
by the explicit formulas 

PI = A I (zi - ^2)-^ + {Z2 - Ci) f + (Ci - C2) f "'f ) , (3.7) 
V q2Pi q^A-ipi q^A-2p^/ 



q^A-^ qt A- 

q2Pi q2A-ipi q;^A-2pi 



={iz2- z,)^ + [z, - C2) + (C2 - Ci) ?r " ) A. (3.8) 



(in) Equations (3.5-3.6) are the discrete Euler-Lagrange equations (1.2) with the Lagrangian function 

£(X, Y) = {Z2 - z,) log(xtxi) + {z, - C2) log(xt A- Vi)+ (3.9) 

(C2 - Ci) log(yjA-Vi) + (Ci - ^2) log(ytA-ixi), (3.10) 

where X = (xijXj), Y = (yi,y2) S Q = x (P''~-^)t. Every vector is an actual r-vector that we 
consider up to resettling. This makes L defined up to an additive constant, but this does not affect the 

Euler-Lagrange equations. 

Proof. Parts (i) follows ininicdiatcly from equations (2.9-2.12) in Theorem 2.2. To establish parts (ii) 

djC 

and (iii), consider the conjugated momentum, P = (6|,62) = ^^(QjQ)- The discrete Euler-Lagrange 
equations (1.2) then split into two groups, 

^\=^{^^Q) = -^{^^Q^) and i32 = -^(Q,Q) = --^(Q,Q). (3.11) 

The first equation becomes 

qjA"'- o^A^^ „t fit A-1 

6l = (^1 - C2)^^ + (C2 - Ci)^f^ = {z, - .2)^ + (.2 - COrrf^, (3.12) 
q^A-ipi q^A-^pi q^^pi q^A-^pi 

which, on one hand, is equation (3.6), and on the other hand, is the equality of two different expression for q| 
in Bf(z), one coming from B^(z)Bf (z) = Bf (z)B^(z), and the other from 'B^{z)'Bf{z) = Bf (^)B^(z). 
The other equation is similar, and that completes the proof. □ 

4. The Isomonodromic Case 

In this section we consider in detail an example of an elementary isomonodromy transformation defined 
in Section 3 of [Kri04]. We show that this transformation can be written in the Lagrangian form and then 
verify that, similarly to the polynomial case considered in [AB06], for rank r = 2 matrices 'L{z) whose divisor 
has 2 simple zeroes and 2 simple poles, this transformation, when written in the spectral coordinates p and 
q, reduces to the difference Painleve equation dPV of the Sakai's hierarchy [SakOl]. Thus, we establish that 
dPV can be written in the Lagrangian form. 

4.1. The Spectral Coordinates. Let L{z) G M^, where V = J2i{zi~Ci)^ and all points are finite, distinct 
and do not differ by an integer. Since the isomonodromy equations (1.11) are invariant w.r.t. the conjugation 
action of the gauge group GLr(C), we can use this action to diagonalize Lg, Lq = diagjpi, . . . ,Pr}, which 
reduces the gauge group to the subgroup C GL^ of the diagonal matrices. Next, consider some asymptotic 
properties of L(^). Namely, let us first introduce the matrix L^o := — reSoo L(^;) ^2: = JZk'^k and put 
ki := ^(L(x))ii = (LQ^L(x))ii. Following [AB06], let us define the type 9 of'L{z) as follows. 

Definition 4.1. The type 9 of the matrix L(2) is the following collection of parameters: 

9{L{z)) = {Zi, . . . , ZnXl, . . . ,Cn; Pi, ■ ■ ■ , Pr, ki, . . . ,kr} . (4.1) 

Prom the multiplicative representation (1.6) we see that these parameters are not independent, since 

n n 

k^ + ... + kr = tr Lo = ^ tr(G,A-i) = ^(z, - Q. (4.2) 

i=l i=l 

For a general choice of parameters, this is the only relation. We denote by Alf the space of matrices of type 
9 and rank r. This space is clearly invariant under the conjugation by non-degenerate diagonal matrices. 
Factoring out this action we obtain the coarse moduli space M^. 

It is worth mentioning that fixing the type 9 of L(z) corresponds to considering a symplectic leaf of the 
canonical foliation of w.r.t. the universal algebro-geometric symplectic form oj of Krichever and Phong 
[KP97, KP98], see the survey [DKP02] for details. 



Lemma 4.1. The dimension of the big cell of the moduli space A4^ is 

diniA^^ = 2(n- l)(r- 1). (4.3) 

Proof. Consider the multiplicative representation of L(z). Each elementary divisor is given by 2(r — 1) 
parameters. Fixing the diagonal elements of Lqo imposes r — 1 (in view of (4.2)) additional conditions, and 
therefore 

dim A^^ = n(2r - 2) - (r - 1) = (2n - l)(r - 1). (4.4) 
Further action by the diagonal matrices reduces the dimension by r — 1, and so 

dimA^f = 2(n-l)(r-l). (4.5) 

□ 

Thus, when n = r = 2, Aif is a complex surface. The spectral coordinate system {q,p) on the space M.2 
is given by the zero q of L(z)i2 and the value of L(z)ii at q normalized in the following way: 

P=- — 7^L(5)ii, where L(g)i2 = 0. (4.6) 
Q- (,2 

The spectral coordinates p and q are essentially the same as in [AB06]. From the integrable systems 
point of view these coordinates are a particular case of the Darboux coordinates of the universal algebro- 
geometric symplectic form cj. For a particular case of the hyperelliptic KdV curves these coordinates were 
first considered by Novikov and Veselov [VN82], the general case was recently established by Krichever in 
[KriOOa, KriOOb]. 

Note that since 

"P(g-C2) 



detL(^) = detLo^^-^^T^^, L(g) 







PiP2(g-Ci) 

p(,q-Z2) 



{z - Zl){z - Z2) 

Lemma 4.2. The additive representation of'L{z) in the pq- coordinates is given by 



Li 



Lo = 



q- zi 

Z2 - Zi 

q- Z2 

Zl - Z2 



P2{q - Z2 + k2) 



P2iq - zi + k2) 



K9-C1) 



Pi{q - Z2 + ki) 



Pip2iq-zi)(q-Ci) 

p{q-Z2) 



[pi{q -zi + ki) - p{q - C2) 1] 



(4.7) 

(4.8) 
(4.9) 



Proof. By choosing a linear normalization in which the first components of the column vectors and the 
second components of the row vectors are equal to 1 , we get 



Hz) 

ing eq 
L(5)i2 = 



Pi 
P2 



[6x 



Z — Z\ 

We then obtain the following equations for ai, 0.2^ bi, 62: 

«! a2 



+ 012 



1 



[62 



Z2 



P = 



q- Zl 

q- Zl 

9- C2 

q- Zl 
9- C2 



+ 



Z2 



L((7)ii = 
Pi + 



Pi 



ai 



q- Zl 



<2 

-(61-62) 



Zl 



-hi 



Oi2 



(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 

Note that we have only three equation for the four unknowns. This is due to the conjugation action by the 
constant diagonal matrices D = diag{7i,72}. Taking this action into account. 



q- Zl 

q- 

ai 



q-Z2 



(Lcx))ll = OL\bi + Q.2b2 = Piki. 



D(aiaibt)D-i = 



71 
— ( 

72 



(4.15) 



we see that we can provisionally put ai = 1. Then 

q- Z2 

0.1 = 

bi = piki - a2&2 = Piki + - — 



Z2 



62 = 62 + (g - C2)P - (a - zi)pi, 



and so 



b2 

bi 



q- Zl 

Z2 - Zl 

q- Zl 



{pi{q - Zl + ki) - p{q - C2)) , 
{pi{q - Z2 + ki) - p{q - C2)) • 



Z2 - Zl 

To find and we use the definition of k2 and detL(g): 

(Loo)22 = Qfia^ + 0:20^ = ^2^2 



detL(g') = det 



9-C2 

q-zi 
* 



P 



P2 



= P1P2 



and so 



Thus, 



q-zi 

(g-Ci)(9-C2) 
{q- z-i){q- Z2y 

P1P2 ^ [q - C,i){q - Zl) 
P q-Z2 







q-z2 



Zl 



P2 



Zl 



P2{q- Zl). 



1 , 2 , , q-^2 2 2 , P1P2 {q-Ci)iq-zi) 

a = p2k2 - a.2a = P2k2 H a = a + 



Zl 



q- Z2 



- p2{q-zi), 



which results in 

,1 



a = 



Zl 



Z2 - Zl 

q- Zl 

Z2 - Zl 



P2{q-Z2 + k2)- 



PlP2(g- Cl) 
P 



P2{q- Zl + k2) 



P\P2 {q - C,i){q - Zl) 



q- Z2 



Conjugating by the diagonal matrix D = diagjq — zi, Z2 — zi} finishes the proof. 



(4.16) 
(4.17) 

(4.18) 

(4.19) 
(4.20) 

(4.21) 

(4.22) 

□ 



4.2. The Additive and the Multiplicative Representations of L(2;). To obtain the multiplicative 
description of L(^) in the spectral coordinates, we need the following Lemma. 



Lemma 4.3. Let 



L{z) = Lo + + = Bf (^)B^(^), 



Zl 



Z2 



where Lq = A^. Then 
(i) the additive representation is given in terms of the multiplicative representation by 

Li = GiB^(^;i), L2 = Bf-(22)G2; 
(a) the multiplicative representation is given in terms of the additive representation by 

- , (^2 - C2) - tr(Lo IL2) 



(zi-Ci)^tr(Lo-%) 
tr(Lo L1L2) 



AG2 



tr(Lo^LiL2) 



LiL 



1-LI2- 



(4.23) 



(4.24) 



(4.25) 



Proof. Part (i) follows immediately from taking the residues at Zi. To establish part (ii), note that from part 

(i) it follows that 

GiA = Li + -^, AG2=L2 + -^, (4.26) 



Z2 - Zl 



Zl - Z2 
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and so what we really need is to establish the identity 



G1G2 = {Z2 - Zi) 



(zi-Ci)-tr(LoiLi) 



L1L2 = {zi - Z2) 



{Z2 - C2) - tr(Lo-iL2) 



L1L2. 



(4.27) 



tr(Lo^LiL2) tr(Lo^LiL2) 

Since Lj and Gi are of rank one, from (i) we see that G1G2 = 7L1L2. The proportionality constant 7 can 
be found using either of the equations (4.26). For example, using the first equation, we see that 



tr(GiG2A-2) (z2-zi)tr(GiA-i-LiA-2) (zi - Ci) - tr(Lo ^Li) 

7 = , ^ ^ = , ^ X = (-22 - 21)- 



tr(Lo^LiL2) 



tr(LoiLiL2) 



tr(Lo'LiL2) 



Corollciry 4.4. The multiplicative representation of'L{z) in the spectral coordinates is given by 



where 



L{z) 



GiA = 



GiA 

z — Zi 



P 



AGo 



P-Pi 

Pi 

Pi -P 



Lq ^ ( Ln 



1 



p2iq - Z2 + k2) - 



-P2{q-k2-Ci) + '^{q-z2) 



AG2 

Z - Z2 



P 

P2P( 



Pl(«-fcl-C2)+'^^ 1 

[Pliq + ki - zi) - p{q - (2) 1 



Proof. Proof is a direct calculation. 



(4.28) 

□ 

(4.29) 

(4.30) 

(4.31) 

□ 



4.3. The Isomonodromic LagreingiEin £tnd dPV. We are now ready to prove the main Theorem of this 
section. 



Theorem 4.5. Consider the special isomonodromic transformation of the form 

Mz) = Bf (z)B^(z) ^ Uz) = (z)B^(z) = B^iz + l)Bf (z) = B^(z + l)L(z)(B^(z))- 
(i) This transformation satisfies the time- dependent Euler-Lagrange equations 

^(Qfe-i, Qfe, fc - 1) + ^(Q/c, Qfe+i, fc) = 
with the Lagrangian function given by 

£{X, Y, t) = {Z2{t) - zi) log(x|xi) + (zi - C2(t)) log(x|A- Vi)+ 



(4.32) 
(4.33) 

(4.34) 
(4.35) 



(C2(t) - Ci) log(yiA-Vi) + (Ci - Z2{t)) log(y|A-ixi), 
where Z2(t) = Z2 — t and C2(i) ~ C2 — t. 
(ii) This transformation defines a birational map ip '■ ^2 ~^ ^2> where Z2 = Z2 — 1, C2 = C2 — 1 o,nd 
all other parameters are unchanged. In the spectral coordinates p, q on M.2 and p, q on M2 
transformation is given by the difference Painleve equation dPV of the Sakai's hierarchy, 

Pl{k2 + Cl-Z2) , P2{k2 + C2-Zi) 

(4.36) 



Z2 + C2 



P- Pi 

{q- Ci){q- zi) 

PP = P\P2 — 



+ 



P- P2 



' {q - Q2){q - Z2) 

Proof. Part (i) is proved in exactly the same way as Theorem 3.1 part (iii). To establish part (ii), it suffices 
to rewrite equation (4.32) in the form (B^(z + l))~"'^L(z) = Bf"(z) and then evaluate it at q. Prom the 
definition of q and p, we get 



Lq ^ I Ln 



AG, 



A-i 



P(9-C2) 





PiP2{q-Ci) 



GiA 

q- zi 



(4.37) 



AG \ / G A 

Thus, second columns of the matrices ^ I Lq — ^3 — ) and ( Lq + 1 are proportional with the 



9- C2 



q- Zl 



proportionality coefficient ^^^1^^ - ^^'^ • Direct calculation then completes the proof. 



p{q - Z2) 



□ 
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